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Transport coefficients associated with the mass flux of an impurity immersed in a granular gas 
under simple shear flow are determined from the inelastic Boltzmann equation. A normal solution 
is obtained via a Chapman-Enskog-like expansion around a local shear flow distribution that retains 
all the hydrodynamic orders in the shear rate. Due to the anisotropy induced by the shear flow, 
tensorial quantities are required to describe the diffusion process instead of the conventional scalar 
coefficients. The mass flux is determined to first order in the deviations of the hydrodynamic 
fields from their values in the reference state. The corresponding transport coefficients are given in 
terms of the solutions of a set of coupled linear integral equations, which are approximately solved 
by considering the leading terms in a Sonine polynomial expansion. The results show that the 
deviation of these generalized coefficients from their elastic forms is in general quite important, even 
for moderate dissipation. 
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I. INTRODUCTION 

The study of transport properties of an impurity or intruder immersed in a granular gas described by the inelastic 
Boltzmann equation is quite an interesting problem. In particular, when the gas is in homogeneous cooling state 
(HCS), the mass flux j for the impurity can be obtained by solving the Boltzmann equation by means of the 
Chapman-Enskog expansion [lj around the local version of the HCS. In the first order of the expansion, the mass 
flux jo is linear in the gradients of mole fraction, pressure and temperature and in the external applied force. The 
corresponding transport coefficients are the diffusion coefficient D, the pressure diffusion coefficient D pi the thermal 
diffusion coefficient Dt and the mobility X- As in the elastic case, these coefficients are the solutions of a set of coupled 
integral equations 0, Q which can be approximately solved by considering the leading terms in a Sonine polynomial 
expansion. Explicit expressions for these transport coefficients have been obtained H[ in terms of the coefficients 
of restitution and the parameters of the mixture (masses and sizes). These analytical results show in general a good 
agreement with those obtained from computer simulations 0, HI, even for strong dissipation. 

However, the above HCS is not accessible experimentally. In order to keep a granular system fluidized, an external 
energy supply is required. In some experimental situations, the gas is driven into rapid flow by the presence of a shear 
field. In this case, a steady state is reached when viscous heating is exactly balanced by the energy dissipation in 
collisions. This unbounded shear flow problem is usually referred to as the simple or uniform shear flow (USF) and 
its study has received a great deal of attention in the past years, especially in the case of monocomponent systems 
[1, Q. Nevertheless, much less is known on transport in sheared granular mixtures [l(| HH, E2) EH- ln particular, 
the understanding of mass transport in granular shear flows is of practical interest, since, for instance, powders must 
frequently be mixed together before any sort of process can begin. Due to the complexity of the general problem 
(binary sheared granular mixture), to gain some insight one can first consider the special case of impurity or tracer 
particles, namely, a binary mixture where the mole fraction of one of the components is negligible. The tracer or 
impurity problem is more amenable to analytical treatment since tracer particles are directly enslaved to the granular 
gas and there are fewer parameters than in a general binary mixture problem. Even in this limit case, the analysis of 
mass transport in a strongly shearing granular gas is an intricate problem basically due to the anisotropy induced in 
the system by the shear flow. For this reason, tensorial quantities {Dij, D p ^j, Dx,ij, Xij} are required to describe the 
mass transport process instead of scalar coefficients {D, D p , Dt, x}- The aim of this paper is to get the above tensors 
in the framework of the inelastic Boltzmann equation. 

Some previous attempts have been carried out earlier in the case of the diffusion tensor D^, especially in the 
self-diffusion problem [141] . However, in all these studies the diffusion was observed in only one direction, usually the 
directionparallel to the velocity gradient. The full diffusion tensor has been obtained in granular gases using kinetic 
theory [Til ], measured in simulations of rapid granular shear [HI ], and even in some experiments of dense, granular 
shear flows in a two-dimensional Couette geometry (l7j . Substantial work on granular diffusivity has been carried out 
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by Hsiau and co-workers [l8(, who have measured self-diffusion coefficients in a variety of granular systems. All these 
results clearly show that diffusion is anisotropic and is significantly affected by inelasticity. 

At a kinetic level, in the tracer limit one can assume that the velocity distribution function /(r, v, t) of the granular 
gas (excess component) obeys a closed (nonlinear) Boltzmann equation, while the velocity distribution function 
/o(r,v,i) of the impurity satisfies a (linear) Boltzmann-Lorentz equation. Both kinetic equations have been solved 
in the (pure) USF problem to get the rheological properties of the mixture in the steady state [13, EH]. Let us 
assume now that the system (granular gas plus impurity) is in a state that deviates from the USF by small spatial 
gradients. In addition, we also assume that the impurity and gas particles are subjected to the action of external 
forces. Since the gas is slightly perturbed from the USF, the Boltzmann equation can be solved by expanding in small 
gradients around the (local) shear flow distribution instead of the (local) HCS. This Chapman-Enskog-like expansion 
has been very recently used [lil , |20| to determine the heat and momentum fluxes to first order in the deviations of the 
hydrodynamic field gradients from their values in the reference shear flow state. Once the state of the excess component 
is well characterized, the goal here is to determine the mass transport of impurity by solving the Boltzmann-Lorentz 
equation by means of a similar perturbation scheme, namely, a gradient expansion around the corresponding shear 
flow distribution which applies for arbitrary values of the shear rate. In the first order of the expansion, the 
tensors Dij, D p _ij, Dx.ij, and Xij are identified from the mass flux jo- As was already pointed out in Refs. (l9l. [20|. 

an important point is that, for general small deviations from the shear flow state, the zeroth-order distribution /q ' is 
not a stationary distribution since the collisional cooling cannot be compensated locally for viscous heating. This fact 
gives rise to new conceptual and practical difficulties not present in the results based on the conventional Chapman- 
Enskog expansion. Due to these difficulties, the results will be restricted here to particular perturbations for which 
steady-state conditions apply and so the (reduced) shear rate (which is the relevant nonequilibrium parameter of the 
problem) is coupled to dissipation. 

The plan of the paper is as follows. In Sec. |TT] we describe the problem we are interested in. The USF state is 
analyzed and the nonzero elements of the pressure tensor for the gas and the impurity are obtained by means of 
Grad's method [H, E EH- Section [□] deals with the perturbation scheme used to solve the Boltzmann-Lorentz 
equation for the impurity to first order in the spatial gradients. The results show that the generalized transport 
coefficients Dij, D p ^j, Dr.ij, and Xij are the solutions of a set of coupled linear integral equations. A Soninc 
polynomial approximation taking the shear flow distribution f{f^ as the weight function is applied in Sec. HVI to solve 
these integral equations and get explicit expressions for these transport coefficients. The details of the calculations 
are displayed along several Appendices. The dependence of some of these transport coefficients on dissipation is 
illustrated in the three-dimensional case, showing again that the influence of inelasticity on mass transport is quite 
significant. Finally, the paper is closed in Sec. |V] with some concluding remarks. 

II. DESCRIPTION OF THE PROBLEM 

As already said, the aim of this work is to analyze mass transport of a dilute granular mixture subjected to uniform 
shear flow (USF) in the tracer limit, i.e., when the mole fraction of one of the species is negligible. In this case, 
the state of the granular gas (excess component) is not affected by the presence of the tracer particles. At a kinetic 
theory level, this implies that the velocity distribution function of the excess component obeys a closed nonlinear 
Boltzmann equation. Moreover, the mole fraction of the tracer particles is so small that their mutual interactions can 
be neglected as compared with their interactions with the particles of the excess component. As a consequence, the 
velocity distribution function of the tracer particles satisfies a linear Boltzmann-Lorentz equation. This is formally 
equivalent to study an impurity or intruder in a dilute granular gas, and this will be the terminology used here. Let 
us start by offering a short review on some basic aspects of the inelastic Boltzmann equation and its solution in the 
USF state. 

We consider a granular gas composed by smooth inelastic disks (d — 2) or spheres (d = 3) of mass m and diameter 
a. The inelasticity of collisions among all pairs is accounted for by a constant coefficient of restitution a (0 < a < 1) 
that only affects the translational degrees of freedom of grains. We also assume that the particles feel the presence 
of an external conservative force F (such as a gravity field). At low density, the time evolution of the one-particle 
velocity distribution function of the gas f(r,v,t) is given by the inelastic Boltzmann equation [221 [23| 

(j t + v • V + ^ • ^ /(r , v , t ) = J[v|/(t), /(*)], (1) 
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where the Boltzmann collision operator J[v|/, /] is 

J[vi\fj} = o d - x I dv 2 J d*e(*-g)(*-g) 

x [a- 2 /(r, vi)/(r, v' 2 ,t) - /(r, v x , t)/(r, v 2 , i)] . (2) 

Here, <r is a unit vector along their line of centers, is the Heaviside step function, and g = vi — v 2 is the relative 
velocity. The primes on the velocities denote the initial values {v'^v^} that lead to {v 1; v 2 } following a binary 
collision: 

vi = Vl - ~ (1 + a- 1 ) (a ■ &)&, V 2 =w 2 + ^(l + a- 1 )(a-g)a (3) 
The first five velocity moments of / define the number density 

n(r,t) = J dv/(r,v,t), (4) 

the flow velocity 

u(r,t) = -±-r f dvv/(r,v,t), (5) 
n(r,t) J 

and the granular temperature 

^ < ^ •^4b)/' ivl ' 2/<^ • v • ^, • (6) 

where V(r, t) = v — u(r, t) is the peculiar velocity. The macroscopic balance equations for density n, momentum mu, 
and energy ^nT follow directly from Eq. ([T]) by multiplying with 1, mv, and ^mv 2 and integrating over v: 

D t n + nV • u = , (7) 
D t u+ (mny 1 (V • P - nF) = , (8) 



AT + A (V • q + PijVjUi) = -(T , (9) 
an 



where D t — dt + u • V is the material time derivative. The microscopic expressions for the pressure tensor P, the heat 
flux q, and the cooling rate £ are given, respectively, by 

P(r,t) = J dvmVV/(r,v,£), (10) 
q(r, t) = J&v ^mV 2 V /(r, v, t), (11) 

C(M) = -^^/dv m yV[ r) v|/(t)]. (12) 

Let us suppose now that an impurity or intruder of mass mo and diameter gq is added to the gas. As said before, 
the presence of the intruder does not perturb the state of the gas, so that its velocity distribution function is still 
determined by the Boltzmann equation {1} . In addition, the macroscopic flow velocity and temperature for the mixture 
composed by the gas plus the impurity are the same as those for the gas, namely they are given by Eqs. ([5]) and 
©, respectively. The velocity distribution function /o(r, v, t) of the impurity satisfies the linear Boltzmann-Lorentz 
equation 

+ v . V + ^ • A) /o(r,v,t) = J [v|/o(*), /(*)], (13) 
at mo av / 
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where the collision operator Jo[ v l/o, /] is now 

•Mvil/o,/] = a 4 ' 1 J dv 2 J d5 0(?-g)(ff -g) 

x [a 2 /o(r,v' 1 )/(r,v^<)-/ (r,v 1 ,i)/(r,v 2 ,t)] . (14) 

Here, Fo denotes an external force acting on impurity, a = [a + <7o)/2 and ao (0 < ao < 1) is the coefficient of 
restitution for impurity-gas collisions. The precollisional velocities are given by 

[l + a Q 1 ) (CT-gi 2 )CT, 



m + mo 



^2 = v a + -5 9 -(l + ao 1 )(5-gi2)^. (15) 
m + mo 

As shown in Ref. [24(, the operator Jo [v|/o, /] is the same as that of an elastic impurity (ao = 1) with an effective 
mass mj = mo + (mo + m)(l — ao)/(l + ao). 
The number density for the impurity is 

no(r,t) = J dv/ (r,v,t). (16) 

The impurity may freely exchange momentum and energy with the particles of the gas and, therefore, these are not 
invariants of the collision operator Jo[v|/o,/]. Only the number density no is conserved, whose continuity equation 
is directly obtained from Eq. (fT5|) 

D t n + n V'UH — = 0, (17) 

m 

where jo is the mass flux for the impurity, relative to the local flow u, 

jo = mo 



i/dvV/o^v,*). (18) 
At a kinetic level, an interesting quantity is the local temperature of the impurity defined as 

This quantity measures the mean kinetic energy of the impurity. As will be shown later, the global temperature T 
and the temperature of the impurity T are in general different, so that the granular energy per particle is not equally 
distributed between both components of the system [l(| [HJ • 

Let us start by describing the state of the system (gas plus impurity) in the (pure) USF. This idealized macroscopic 
state is characterized by the absence of external forces (F = F = 0), constant densities n and no, a uniform 
temperature, and a linear velocity profile u = a • r, where the elements of the tensor a are = aSi X Sj y , a being the 
constant shear rate. This linear velocity profile assumes no boundary layer near the walls and is generated by the 
Lee-Edwards boundary conditions [25j , which are simply periodic boundary conditions in the local Lagrangian frame 
moving with the flow velocity. For elastic gases, the temperature grows in time due to viscous heating and so a steady 
state is not possible unless an external (artificial) force is introduced [26|. However, for inelastic gases, the temperature 
changes in time due to the competition between two (opposite) mechanisms: on the one hand, viscous (shear) heating 
and, on the other hand, energy dissipation in collisions. A steady state is achieved when both mechanisms cancel each 
other and the fluid autonomously seeks the temperature at which the above balance occurs. Under these conditions, 
in the steady state the balance equation © becomes 

aP xy = -~(p, (20) 



where p = nT is the hydrostatic pressure. The balance equation (|20[) shows the intrinsic connection between the 
shear field and dissipation in the system. As a consequence, the shear flow state associated with (|2T))) is inherently 
beyond the scope of the Navier-Stokes or Newtonian hydrodynamic equations [2l[ since the collisional cooling (which 
is fixed by the mechanical properties of the particles) sets the strength of the velocity gradient in the steady state. 
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Furthermore, note that for given values of the shear rate a and the coefficient of restitution a, the relation (|20|) 
gives the temperature T in the steady state as a unique function of the density n. In fact, the reduced shear rate 
a* = a/i/(n, T) is only a function of the coefficient of restitution a in the steady state. Here, v(n, T) is a characteristic 
collision frequency given by 

^ (d ~ 1)/2 8 

The relevant transport properties of the system in the steady USF are related to the pressure tensor P since 
q = jo = 0. Furthermore, one can introduce the partial pressure tensor P defined as 

Po = y dvm VV/ (V). (22) 

Explicit expressions for the nonzero elements of P and Po have been recently obtained from the Boltzmann equation 
by means of Grad's method [HI. A brief summary of these results is given in Appendix lAl The accuracy of this Grad's 
solution has been confirmed by comparison with Monte Carlo simulations of the Boltzmann equation [T3, H3, HH , 
even for strong dissipation. 

III. CHAPMAN-ENSKOG-LIKE EXPANSION: TRANSPORT PROPERTIES AROUND USF 

The main goal of this paper is to determine the mass flux for the impurity in the presence of USF. In that 
case, let us assume that the USF state is disturbed by small spatial perturbations. The response of the system to 
those perturbations gives rise to contributions to the mass flux, which can be characterized by generalized transport 
coefficients. This Section is devoted to the evaluation of these coefficients. 

In order to analyze this problem we have to start from the Boltzmann-Lorentz equation (| 1 3j> with a general time 
and space dependence. Let u s = a • r be the flow velocity of the undisturbed USF state. In the disturbed state, however 
the true velocity u is in general different from u s since u = u s + Su, Su being a small perturbation to u s . As a 
consequence, the true peculiar velocity is now c = v u = V Su, where V = v — u s . In the Lagrangian frame 
moving with u s , the Boltzmann-Lorentz equation can be written as 

Ft fo ~ aVy W x fo + (v + Us) ' v/o + S ' ^v fo = Jo[V|/o ' /L (23) 

where here the derivative V/o is taken at constant V. The macroscopic balance equations associated with this 
disturbed USF state follows from the general equations ((8]), (|9|), and (fTTj) when one takes into account that u = u s +5u. 
The result is 

8 t n + u s • Vn = -V • (n Su) - (24) 

m 

d t Su + a ■ Su + (u s + Su) ■ V<5u = -(mn)" 1 (V • P - nF) , (25) 

^nd t T + ^n(u s + Su) ■ VT + aP xy + V ■ q + P : VSu = ~pC, (26) 

where the pressure tensor P, the heat flux q, the cooling rate C, and the mass flux jo are defined by Eqs. (|10[) . pip . 
(fl"2|) . and (fl"8|) . respectively, with the replacement V — > c. 

We assume that the deviations from the USF state are small. This means that the spatial gradients of the hy- 
drodynamic fields A(r,t) are small. As noted in Ref. [2], there is some flexibility in the case of granular mixtures in 
the representation of the heat and mass fluxes since they can be defined in a variety of equivalent ways depending 
on the choice of hydrodynamic gradients used. In fact, some care is required in comparing transport coefficients in 
different representations using different independent gradients for the driving forces. Here, as in previous works 0,3, 
I take the concentration of impurity £0 — n o/ n y the pressure p, the temperature T, and the local flow velocity Su as 
hydrodynamic fields, i.e., 

A(r,t) = {x (r,t),p(v,t),T(r,t),Su(r,t)}. (27) 
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Since the system is strongly sheared, a solution to the Boltzmann-Lorentz equation ([23)) can be obtained by means of 
a generalization of the conventional Chapman-Enskog method [l[ where the velocity distribution function is expanded 
about a local shear flow reference state in terms of the small spatial gradients of the hydrodynamic fields relative to 
those of USF. This is the main new ingredient of the expansion. This type of Chapman-Enskog- like expansion has 
been considered in the case of elastic gases to get the set of shear-rate dependent transport coefficients [26|, |29[ in a 
thermostatted shear flow problem and it has also been recently considered [HI, H3] for inelastic gases. 
In the context of the Chapman-Enskog method [l[ , we look for a normal solution of the form 

/o(r,V,t)=/ (A(r,t),V). (28) 

This special solution expresses the fact that the space dependence of the reference shear flow is completely absorbed 
in the relative velocity V and all other space and time dependence occurs entirely through a functional dependence 
on the fields A(r, t). Moreover, in the presence of external forces it is necessary to characterize the magnitude of these 
forces relative to the gradients as well. Here, it is assumed that the magnitude of the external forces F and Fo is 
of first order in perturbation expansion. The functional dependence (|28l) can be made local by an expansion of the 
distribution function in powers of the hydrodynamic gradients: 

f (A(r, t), V) = /< 0) (V) + / (1) (V) + • • • , (29) 

where the reference zeroth-order distribution function corresponds to the USF distribution function but taking into 
account the local dependence of the concentration, pressure and temperature and the change V — > V — <5u(r, t) = c. 
The successive approximations are of order k in the strength of the external forces as well as in the gradients 
of Xq, Vi T, and Su but retain all the orders in the shear rate a. Here, only the first order approximation will be 
considered. 

When the expansion (|29|) is substituted into the definitions (fT0|) . (TTTj) . (Tl2|) . and (fP8|) . one gets the corresponding 
expansions for the fluxes and the cooling rate: 

P = P (o) + P (i) + . . . i q _ q (o) + q (i) + . . . ) ( 30a ) 

C = C (0) + C (1) + • ■ • , jo=j^ 0) +j^ 1) + --- - (30b) 

Finally, as in the usual Chapman-Enskog method, the time derivative is also expanded as 

d t = df ] + d { t 1] + 9 t (2) + • • • , (3i) 

where the action of each operator is obtained from the hydrodynamic equations (f2"4"l) - ([2"T)) . These results provide 
the basis for generating the Chapman-Enskog solution to the inelastic Boltzmann-Lorentz equation (|23|) . 

A. Zeroth-order approximation 

Substituting the expansions (f29|) -l[3"l" j) into Eq. (|23l) . the kinetic equation for is given by 

Oi 0) C-aV y ^-C = J [V\Cj% (32) 

where use has been made of the fact that Fo is of first order in gradients. To lowest order in the expansion the 
conservation laws give 

9 t (0) * = 0, r- 1 9 t (0) T = p- 1 9(%=-A aP (o)_ c (o) 7 (33) 



d^'Su. + a.jSuj =0. (34) 

As shown in Refs. .19, 20] , for given values of a and a, the steady state condition (|20|) establishes a mapping between 
the pressure p and temperature T so that every pressure corresponds to one and only one temperature. Since the 
pressure p(r, t) and temperature T(r,t) are specified separately in the local USF state, the viscous heating term 

a|-Pi^| only partially compensates for the collisional cooling and so, the pressure and temperature depend on time. 
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This implies that the zeroth-order distributions for the gas and the impurity depend both on time through 
their dependence on the pressure and temperature. Consequently, in general the reduced shear rate a* — a/v(p,T) 
depends on space and time so that, a* and a must be considered as independent parameters for general infinitesimal 
perturbations around the USF state. This fact gives rise to conceptual and practical difficulties not present in the 
case of elastic collisions ;29j . 

Since /q is a normal solution, the time derivative in Eq. (|32j) can be represented more usefully as 

where in the last step we have taken into account that depends on Su only through the peculiar velocity c. 
Substituting (|35f into |32|) yields the following kinetic equation for : 

To solve Eq. (f36|) one needs to know the dependence of the momentum flux Pxy on the pressure p and temperature 
T. A detailed study of this problem has been carried out in Ref. [2l|. The first nontrivial velocity moment of the 

distribution corresponds to the partial pressure tensor P °^ defined as 

P 0) =m J dccc/ (0) (c). (37) 



The temperature of the impurity T (t) can be determined from the trace of P °' . From Eq. (|36|) , one gets 



(-?-<<> + M (pg- + t| ) P<°> + a,,P<°> + % ,P<°> = B,„ (38) 



where <2y = aSi X Sj V and 



B« - m / dc c iCj J [/ (0) , / W]. (39) 

A good estimate of the collisional moment By can be made by considering Grads's approximation [l5l |. In this case, 
Bij is given by Eq. (|B6j) . The steady state solution of Eq. (|38|) is also displayed in Appendix [A] However, in general 



the equations ([38]) must be solved numerically to get the dependence of the zeroth-order pressure tensor P^}- (p, T) on 

(0) a „A p(°) 



pressure and temperature. The behavior of the pressure tensors and P / near the (steady) USF state is studied 
in Appendix [El In what follows, P^ and Pq^j will be considered as known functions of p and T. 



B. First-order approximation 



The analysis to first order in the gradients is similar to the one carried out in Ref. [20j for the one-component case. 

•(i 
o 



Some details of this derivation are given in Appendix [Cj The distribution function /i 1 ^ is of the form 



= Aq ■ Vx + B ■ Vp + C • VT + V : V5u + £ -F, (40) 

where 

^ = F -^F. (41) 

777, 
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The vectors {Aq, Bo, Co, £q}, and the tensor T)q are functions of the true peculiar velocity c. They are the solutions 
of the following linear integral equations: 



^aPj°> + C (0) ) ( P d P + Td T ) Ao - acy^-Ao - J [Ao, / (0) ] = A 



(42) 



aPg? + C (0) ) (pd P + Td T ) Bo - {^d p P^ + 2^ + ac y ^j B 



-Jo{BoJ {0) ]=B ~ 



^(l-^)Pj°) 



C + J [/o {0) ,B], 



(43) 



dp 



aPf$+C {0) )(pd P + Td T )Co 



2 a 



<1 ^ + Td T )P^ + \^ 



+ac. 



dc x 
(o) 



Co-Jo[CoJ {0) ] = Co+{^d T P^- P - 



< (0) 
2T 



Bo 



+Mfr,c\, 



(44) 



^aP$ + C (0) ) {pdp + Td T ) V 0M - acy^-Vo^ - aS ey V , xj 
-JoPomJW] = Do 4j + C«,« (pd P + Td T ) /C°) + Jol/o^,^], 



(45) 



^aP$ + C (0) ) (pdp + Td T ) £o - acy^£o - Jo[£o, / (0) ] = E . 



(46) 



Here, A (c), B (c), C (c), D (c), and E (c) are defined by Eqs. (jCSJ), (|CT0f . ([UT!]) . (jCT2| . and (jCT3)l . respectively. 
In addition, upon writing Eqs. (|4*2")) - (j4*6"|) , use has been made of the explicit form of /W. It has been derived in Ref. 
12011 and reads 



= B ■ Vp + C ■ VP + T> : V5u, 



(47) 



where the coefficients B, C, and T> are functions of the peculiar velocity c and the hydrodynamic fields. In particular, 
the first-order contribution to the cooling rate, Cuiji is given by 



~ J dcmc 2 (j[/ (0) ,Ai] + J[%,/ (0) ] 



The first-order contribution to the mass flux Jq of the impurity is defined as 

Jo 1} = m J dcc/ (1) (c). 
Use of Eq. (|4*0~| into Eq. (jl9")l gives the expression 



.(i) _ dx m dp mn dT 



where 



d rj T p ' lJ dr 3 T 



D i: j = - J dec,; A j(c), 



(48) 



(49) 



(50) 



(51) 



Tm 



m 



J dcci B ,j(c), 



(52) 
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D T ,ij = 



Tmi 



■ J dccj C ,j(c), 



(53) 



Xij = m a / dec; fo,i(c). 



(54) 



Upon writing Eqs. (f5"Tj) -(f54 f use has been made of the symmetry properties of *Aq, Bo, Co, and So. In general, the 
set of generalized transport coefficients -Dy, D p ^j, Dx,ij, and Xij are nonlinear functions of the shear rate and the 
coefficients of restitution a and ao ■ It is apparent that the anisotropy induced by the presence of shear flow gives rise 
to new transport coefficients for the mass flux, reflecting broken symmetry. According to Eq. (f50|) . the mass flux of the 



impurity is expressed in terms of a diffusion tensor D^j , a pressure diffusion tensor D p ^j , a thermal diffusion tensor 
Dx : ij, and a mobility tensor Xij- Note that in the particular case of the gravitational force F = mg and F = m g, 

Consequently, 



where g is the gravity acceleration. In this case, the combined force T defined in Eq. (|4ip vanishes, 
the external force T does not occur in Eq. ([40)1 when the system is only subjected to a gravity field. 



C. Steady state conditions 



The evaluation of the above transport coefficients requires to know the complete dependence of pjp and Pq ^- on the 
pressure p and the temperature T. This involves the corresponding numerical integrations of the differential equations 



obeying the pressure tensors Pjp and Pq U ^ [see Eq. (j3"8")) for Pq^]- Needless to say, this is quite an intricate problem. 
However, some simplifications occur if attention is restricted to linear deviations from the USF steady state described 
in Sec. HIl In particular, since the contributions to the mass flux (|50[) are already of first order in the deviations from 
the steady state, one only needs to know the transport coefficients to zero order in the deviations. This means that 



>(0) 



>(0)l 



9 t w T = df ] p = and so the term (2/dp)aP^y > + C (0) = in the left hand side of Eqs. gj)-®- The differential 



,(0) 



equations for the transport coefficients thus become simple coupled algebraic equations. 

on the pressure p and temperature T occurs explicitly and through its dependence on the 
oc \pT jp. Consequently, 



The dependence of P\ 
reduced shear rate a 



pd p P^ ] = pd p pP: 3 {a*) = p ( 1 - a*-^ ) P*(a*) 







(55) 



Td T P^ = Td TP P*(a*) = lpa*^P*(a*) 



(56) 



where P*j = P^ /p and a* — a/ v. The dependence of P*j on a* near the steady state was determined in the one- 
component problem [20| so that all the terms appearing in the integral equations are explicitly known in the steady 
state. Under the above conditions, Eqs. (|42 p -(|46 |) become 



d 

- aCy—Ao - Jo[Ao, / (0) ] = A , 



(57) 



2a 
~~d 

2aT 

—a d a ,k> xy 



I - a*d a *)P* y + 2C (0) 



ac. 



dc, r 



Bo- Jo [B ,f 



(0)1 



Bn 



Co + Jo[/ (0) ,S] 



(58) 



^(l + ia*0 a .)i* 



d 

p j a 
"f 



2 C +aCy dc x 

p(0) 



Co- Jo[C ,/ (0) ] =C 



d a*d a «P: y - l—)B + Jo[f^,C], 



(59) 



ac y ^—Voi 3 - a6 ey V 0tXj - J [D 0M , /(°)] = D 0M + ( uM { P d p + Td T ) / (0) + Jo [/<T , Vtj] , 

UCx 



(60) 
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- ac y — £ Q - Jo[£o, / (0) ] = E . (61) 

It must be recalled that in Eqs. (f57]) - (|6ip all the quantities are evaluated in the steady state, namely, P*j and a* 
are given by Eqs. (|A3 |> — (|A5 p and (|A6|> . respectively, while dP^/da* is given by Eqs. (|BT|) - (jB3)) . Henceforth, the 
calculations will be restricted to the particular condition (|2T))) . 



IV. MASS TRANSPORT OF THE IMPURITY 



This Section is devoted to the determination of the generalized transport coefficients associated with the mass 
transport of the impurity. In order to get explicit expressions for these coefficients, one has to know the quantities 
Aq, B , C , and £ which verify the coupled integral equations ([ST)) . 1(55]). ((55)) . and (|6"Tj). respectively. To determine 
the explicit dependence of these quantities on the coefficients of restitution a and ao, one needs to make use of certain 
approximations. The standard approach is to consider the leading term in a Sonine polynomial expansion. In a 

previous work on diffusion in shear flow [l5T |. the isotropic part of in this expansion was assumed for simplicity 
to be a Maxwellian /o,m(c). However, given that the system is strongly sheared, it is reasonable to expect that the 

isotropic part of /q is mainly governed by the shear flow distribution /q rather than by the Maxwellian distribution. 
For this reason, here we keep the usual structure of the standard Sonine approximation, except that the Maxwellian 
weight function /o,m is replaced by A similar type of modified Sonine approximation has been recently considered 
[30] | to estimate the Navier-Stokes transport coefficients of a single granular gas. 

According to the above arguments, in the case of the mass flux jg , a good estimate of {Aq, Bq, Co, £q} is given 
by the first Sonine approximation: 

A),i ~ c 3 D 3i /o ( c ); &o,i -> -Cj D p tji flf\c), (62a) 



Co,,: - <■„!),, ., f$» (c), £ 0l , - / W (c), (62b) 

where fo°\c) is the solution of Eq. (|36|) . The relationship between the tensors {-Dy, Dp,ij, Dr.ij,Xij} an d their 
corresponding counterparts {D^ , D p t j , P>T,tj > Xij} i s consistently obtained from Eqs. (|51[) - (|54[) . A simple calculation 
yields 

Dij = moQikDkj, Dp : ij = 7fiQikDp,kji (63) 



— 111/111 

fr,ij = —QikD T ,kj, Xij = QikXkj, (64) 

where Q = Pq 1 . When one takes /o,m instead of /q "* in the Sonine expansion, then = (n-o7o) _1 ^ij and one 
recovers previous results (15j . Consistently, the quantities S and C corresponding to the distribution /W of the 
granular gas [see Eq. (|47[) ] must be similarly approximated. However, as shown in Ref. [20(, these quantities vanish 
in the first Sonine approximation and so there is no contribution to the mass flux coming from the terms of the 
form J [f { ° ] ,&] and J [f^ ] ,C]. Substitution of the expressions ([52"a]) and ([62bf into Eqs. ([57 ]) . ([55 ]) . ([55 ]) . and ([61 ]) 
gives a closed set of integral equations for D^, D v .ij, Dt.ij, and Xij- Multiplication of these equations by moc, and 
integration over c yields 

P 

(a lk + Q ik ) D kj = — P ij , (65) 



^(l-a*d a ,)P: y + 2^ 



(") 



D 



(an- + fijfe) D Pt kj — 

m \ m J J 

2 -l a *d a *P* xy -^)D T , 3 , 



(66) 
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FIG. 1: Plot of the reduced elements (a) D* x , (b) D% z , (c) D* y , (d) D yx , and (e) D* y as functions of the (common) coefficient 
of restitution a = ao for a three-dimensional system in the case a — a and mo/m = 0.5. 
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FIG. 2: Plot of the reduced elements (a) Dp ZZ , (b) D'p tVy , (c) Dp >xx , (d) Dp >yx , and (e) D* v ^ xy as functions of the (common) 
coefficient of restitution a = «o for a three-dimensional system in the case a = a and mo/m = 0.5. 



(a^ + Qif.) Dr,kj — 

I1^ a * d . f^P* -P* 
2 m ° V m y 



5 a«0..j» 2 



5^) Am 



(ajfc + f^fc) Xkj = noSi 



In the above equations, P * y = P^/xop and = n.QToAikQkj, where 



A, 



«-oT 



dec, Jofo/o^/^]. 



This quantity can be evaluated by using standard integration techniques with the result 

V2 (a^ 



A,- 



Ad 



ufi(l + a ) [(1 + 0)0]" 1/2 {{d + 2)(1 + 0)<^ 



+0 (P* - + \d + 3 + (d + 2)9} (7- 1 i > *,« " %) } 



(67) 
(68) 

(69) 



(70) 



where z/ is defined by Eq. (|2ip . /i = m/(m + mo), 7 = Tq/T is the temperature ratio and = moT/mTo is the 
mean-square velocity of the gas particles relative to that of the impurity. Upon deriving Eq. (fT0"|) , use has been made 

of the Sonine approximations (|A1|) and (|A8|) for the distributions /^°^ and /q "*, respectively, and the nonlinear term 
proportional to (P*- — Sij) (7 _1 Po* i y — has been neglected. 
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FIG. 3: Plot of the reduced elements (a) xy , (b) D^ xx , (c) D^ yx , (d) D^ yy , and (e) D^ zz as functions of the (common) 
coefficient of restitution a = ao for a three-dimensional system in the case a = ~a and mo/m — 0.5. 
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FIG. 4: Plot of the reduced elements (a) x*z> (b) Xyyi ( c ) XaUj (d) Xyx j an d (e) XaSy as functions of the (common) coefficient 
of restitution a = ao for a three-dimensional system in the case a — a and mo/m = 0.5. 



The coefficients Dij and Xij decouple from the other ones and hence can be obtained more easily. By using matrix 
notation, they are given by 

D = ^(a + fi)- 1 .P*, (71) 
m 



X = n (a + n) (72) 

The remaining coefficients D p ^j and Dx,ij are coupled and they obey the set of simple algebraic equations (|66p and 
(EH). 

In the elastic limit (a = ao = 1, which implies a* = in the steady state conditions), T = Tq, P,*j = P *y = 8ij, and 
Slij = floSij, so that Eqs. (frj5|) - (frj5)) have the solutions = D^5ij, D p ^j = D p o5ij, Dx,ij = 0, and Xij = Xo^iji where 
Do, D p o and %o ar e the conventional Navier-Stokes transport coefficients for ordinary gases hfl. Their expressions are 



Ai = D 



where 



Tx / m \ n 

I 1 )> Xo = tt-, (73) 

o v m / sZq 



4 7r (d-i)/2 _ d _ i / 2Tm 
^ 0= dT(W nCT " imoim + moY (?4) 

To illustrate the dependence of the tensors Ty = {Dij,D Pi ij,DT,ij,Xij} on dissipation, let us consider a three- 
dimensional system. In this case, according to Eqs. ([65|) -([68 |l . T xz — T zx = T yz = T zy = 0, in agreement with the 
symmetry of the problem. As a consequence, there are five relevant elements: the three diagonal (T xx , T yy , and T zz ) 
and two off-diagonal elements (T xy , and T yx ). The integral equations (|65p - (|68[) also show that T xx ^ T yy ^ T zz and 
T xy ^ T yx . In Figs. HHH we plot the relevant reduced elements of tensors D*p D*^, Dj,^ and \*j as functions of 
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FIG. 5: Plot of the ratios (a) D% x /xxxj (b) -D* z /x*z, and (c) Dy V /Xy V as functions of the (common) coefficient of restitution 
a = «o for a three-dimensional system in the case a = a and rao/m = 0.5. 



the (common) coefficient of restitution a — ag when a = a and mo/m = 0.5. Here, the tensors have been reduced 
with respect to their values in the elastic case, namely, D*j = Dij/D , D* ^ = D Pt ij/D p o, and x*j = Xij/XOi while 
Dj-,ij = D Tii j /{xqT /mu). We observe that in general the influence of dissipation on the transport coefficients is 
quite significant. This means that the deviation of the elements T^- from their functional forms for elastic collisions 
is important for moderate dissipation. It is also apparent that the anisotropy of the system, as measured by the 
differences \T XX — T yy \ and \T yy — T zz \, grows with the inelasticity. This anisotropy is much more important in the 
plane of shear flow (\T XX — T yy \) than in the plane perpendicular to the flow velocity (\T yy — T zz \). 

As expected, the usual Einstein relation between the diffusion and mobility coefficients for ordinary fluids [3X| ] is 
no longer valid in this noneq uilibrium situation for granular gases. A similar conclusion has been found when the 
gas is under HCS [13, HH, |34|. There are basically two independent reasons for this violation: the occurrence of 
different kinetic temperatures between the impurity and gas particles and the inherent non-Newtonian properties of 
the reference state. The deviation of the ratio D*^/x*j from unity is a measure of the violation of Einstein's relation 
in the USF state. This ratio is plotted in Fig. [5] for the diagonal elements in the case a = a and mo/m = 0.5. As Fig. 
[5] clearly shows, the ratio D*j/x*j ^ 1 and quickly decays in this situation (when the impurity is lighter than the gas 
particles) for the yy and zz elements, while the opposite happens for the xx elements. The violation of the Einstein 
relation obtained here contrasts with some numerical experiments performed by Makse and Kurchan [35j who applied 
uniform shear to measure diffusivity and mobility in bidisperse mixtures. Their results show that Einstein's relation is 
verified so that a temperature can be defined in the system by analogy with ordinary fluids. However, this conclusion 
also disagrees with the lack of equipartition in a granular shewed gas when different particles are present 0. 



V. CONCLUDING REMARKS 



In this paper, mass transport of an impurity or intruder immersed in a strongly sheared granular gas at low-density 
has been analyzed. We have been interested in a situation where weak spatial gradients of concentration, temperature 
and pressure coexist with a strong shear rate. In addition, we have also assumed that the system (gas plus impurity) 
feels the action of external forces which are considered to be at least of first order in the spatial gradients. Under these 
conditions, the resulting mass transport process is anisotropic and thus it cannot be described by scalar transport 
coefficients. Instead, it must be described by shear-rate dependent tensorial quantities whose explicit determination 
has been the main objective of this paper. 

In the tracer limit, the inelastic Boltzmann and Boltzmann-Lorentz kinetic equations describe the state of the gas 
and the impurity, respectively. Since the state of the gas slightly deviates from the USF by small spatial gradients, a 
generalized Chapman- Enskog method has been recently proposed [l9l |2C| to analyze transport around nonequilibrium 
states. Here, a similar perturbative scheme has been used to solve the Boltzmann-Lorentz equation to first order in 
the deviations of the hydrodynamic field gradients from their values in the reference shear flow state /q "* . As noted 
in previous works (l9l [2fj| , the zeroth-order distribution is not in general stationary and only in very special 
conditions has a simple relation with the (steady) USF distribution. Since we are mainly interested in determining 

mass transport of the impurity jg in the USF state, for practical purposes the results have been specialized to the 
steady state conditions, namely, when the hydrodynamic variables satisfy the relation (p?0)) . This implies that the 
reduced shear rate a* is coupled with the coefficient of restitution a, so that the latter is the relevant parameter of 
the problem. Under these conditions, Jq 1 '' is given by Eq. (|50[) where the corresponding set of generalized transport 
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FIG. 6: Plot of the ratios (a) D p ) '(^Dp^k), (b) _D/(^Dfcfc), and (c) x/(^Xfcfc) as functions of the (common) coefficient of 
restitution a = ao for a three-dimensional system in the case a = a and rao/m = 0.5. 

coefficients {Dij,D p ^j, Dx,ij,Xij} are the solutions of the linear integral equations (|5"T)) - (|6"l"j) . To get explicit results, 
a first Sonine polynomial approximation has been considered to estimate some collisional integrals. The reliability of 
this approach has been assessed in the USF problem, where it has been shown to agree very well with Monte Carlo 
simulations [HI, [28[ . 

The results show that the coefficients {Dij, D p ^j, DT^j,Xij} present a complex dependence on the coefficients of 
restitution a and ao and on the masses and sizes of the system. This is clearly illustrated in Figs.[TH4j The deviations 
of {Dij , D Pii j , D Tii j , Xij} from their elastic counterparts are basically due to three different reasons. First, the presence 
of shear flow modifies the collision frequency of the elastic diffusion problem £!o [defined by Eq. (|T4"|) ] by the tensorial 
term + f2jj. Second, given that in general the impurity and gas particles are mechanically different, the reference 

shear flow state of the impurity is completely different from that of the gas particles /(°). This effect gives rise to 
terms proportional to (rno/m)/ 3 * — Pq y. Third, there is a coupling between the coefficients D p ^j and Dx,ij due to 
the inherent non-Newtonian features of the USF state of the gas. Each one of the three effects is a different reflection 
of the dissipation present in the system. 

Most of the works 36] on granular mixtures have been based on the CE expansion around an elastic (local) 
equilibrium state up to the Navier-Stokes order, and therefore they are limited to nearly-elastic systems in the 
USF. A more recent CE expansion 0, H3] around the (local) HCS takes into account energy nonequipartition and 
provides expressions for the Navier-Stokes transport coefficients of the mixture without any restriction on the level 
of inelasticity. In particular, the mass transport of impurity is characterized by the single scalar coefficients D, D p , 
Dt, and x- Although the base state considered in Refs. and [37J is different from the one chosen here, it is still 
worthwhile to carry out some comparison between both descriptions. To that end, let us define the scalars ^Dkk! 
2-D p ,kk, and ^Xfcfc- These coefficients can be understood as the generalized diffusion coefficient, pressure diffusion 
coefficient, and mobility coefficient in a strongly sheared mixture. In Fig. [6J the ratios D/ (^Dkk), Dp/^Dp^k), and 
xl (^Xfcfc) are plotted versus the (common) coefficient of restitution a = a.$ for a three-dimensional system in the case 
cr = a and mo/m — 0.5. Here, D, D p and x refer to the coefficients obtained in Refs. [2| and [32j. As expected, 
we observe that in general the transport coefficients of the perturbed USF state differ from the usual Navier-Stokes 
coefficients as the collisions become more inelastic. These discrepancies are specially significant in the cases of the 
diffusion and the pressure diffusion coefficients. 

It is apparent that the results presented here are relevant to make a comparison with numerical simulations. In the 
self-diffusion problem (when the impurity and gas particles are mechanically equivalent), previous results obtained for 
the diffusion tensor [l5| have shown good qualitative agreement with molecular dynamics simulations [l6|. Beyond 
this particular case, to my knowledge no previous studies on the tensors {Dij, D p ^j, Dx,ij, Xij} have been performed. 
I hope that this paper stimulates the performance of such computer studies to check the relevance of kinetic theory 
to describe mass transport under shear flow. 

A possible application of the results reported in this paper is to study segregation induced by a thermal gradient. 
Thermal diffusion is caused by the relative motion of the components of a mixture due to the presence of a temperature 
gradient. Due to this motion, concentration gradients subsequently appear in the mixture producing diffusion that 
tends to oppose those gradients. In the steady state, j$ = and the thermal diffusion factor provides a segregation 
criterion. Recent kinetic theory results [38[ based on the Navier-Stokes transport coefficients have been able to explain 
some experimental and/or molecular dynamics segregation results obtained in agitated granular mixtures at large 
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shaking amplitudes [39] . Another possible direction of study is the extension of the present approach to mixtures with 
finite composition. Given the mathematical difficulties associated with the description of multicomponent systems in 
far from equilibrium situations, one could perhaps use a kinetic model of the Boltzmann equation. Once the transport 
coefficients of the mixture are known, a linear stability analysis of the hydrodynamic equations could be carried 
out to identify the conditions for instabilities at long wavelengths [20]. Finally, it must be noted that the results 
reported here have been made in the context of a very simple collision model where the coefficients of restitution 
a and ao are constant. However, experiments and simulations show that the coefficients of restitution depend in 
general on the relative velocity of colliding particles • Recent results [U [42| derived for these viscoelastic models 
in the case of the Navier-Stokes transport coefficients show qualitative differences from the ones obtained with the 
simplifying assumption of a constant coefficient of restitution. For this reason, it would be interesting to extend the 
present description to this kind of models (where collisions are described by an impact velocity dependent coefficient 
of restitution) in order to check whether the behavior predicted here for the generalized transport coefficients also 
occurs (at least from a qualitative level) in the case of granular gases of viscoelastic particles. As in the Navier-Stokes 
description [HI], explicit results for this more realistic model could be obtained only for small enough dissipation. 
This contrasts with the results reported here since they are not restricted to any level of inelasticity. Work along the 
above lines will be done in the near future. 
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APPENDIX A: RHEOLOGICAL PROPERTIES IN THE STEADY USF 

The explicit expressions for the pressure tensors P*j = Pij/p and Pq ^ = Po^j/xop of the gas and the impurity in 
the steady USF are provided in this Appendix. To get the explicit expressions of the elements of P*j, one takes the 
following Sonine approximation for /: 



/(V) - / M (V) 



2T V nT 



v k v t - -y 2 5, 



(Al) 



where 



}M{Y) = n {-£r) exp v 27 



(A2) 



By using this approximation, the nonzero elements of P*j are given by [2T| 



p* _ p* 

yy zz 



d+l + (d-l)a 
2d + 3 -3a 



(A3) 



d+l + (d-l)q 
xy (l + a)(2d + 3-3a) 2 



(A4) 



P* xx = d-{d-l)P* 



yy 



(A5) 



The relationship between the reduced shear rate a* = a/v [where v is defined by Eq. ([2"Tjl ] and the coefficient of 
restitution a is 



,2 _ d+ 2 (1 + a) (2d + 3 - 3a) 2 (l - a 2 ) 
a ~ ~32d d+ 1 + (d- l)a 

Moreover, the (reduced) cooling rate £* = C*" "* l v is 



(A6) 



(A7) 
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In the case of Pg i7 -, one considers the leading Sonine approximation for given by 



/o(V) -> /o,m(V) 



mo 
2T 



P 



V k Vt 



-y 2 5, 



(A8) 



where 



/o,m(V) = n 



"'0 
2ttT 



d/2 



cxp 



(A9) 



Another different approach based on an anisotropic Gaussian distribution has been considered by Lutsko [13J. This 
latter approximation has the additional advantage of being positive definite. However, the algebraic equations defining 
the pressure tensor Pg.ij cannot be solved explicitly since it requires to numerically solve some collision integrals. Here, 
for the sake of simplicity, we have preferred to estimate the rheological properties of the system by means of Grad's 
solution (|A8I) . This allows us to get explicit expressions for f *ir ^ nas been shown, by comparison to DSMC 
simulations [12L IT3 |. that the results derived from both approaches compare very well with computer simulations, 
even for strong dissipation. Using (|A8| . the nonzero elements of P * ^ can be written as [lj| 



P 



= P< 



0,zz 



P 



F + HR 



0,dd 



[I [I 



G 



(A10) 



p* _ 

0,xy 



p* _ up* 
^O.yy n fx 

G 



(All) 



d 7 -(d-l)P, 



0,yy> 



(A12) 



where 7 = To /T is the temperature ratio and 

VlfaV- 1 (1 + 



1/2 



(I +"<■) l + |(d-l)(l + fl)(l + a ) 



(A13) 



^ _ V2 fa 



d-1 



1/2 



(1 + Qfo) 



x {2[(d + 2)9 + d + 3] - 3//(l + 0)(1 + a )} , 



(A14) 



4d 



d-1 



/*0 



0(1 



1/2 



(l + oo)[3/i(l + 0)(l + a o )-2]. 



(A15) 



Here, = mo/mj, [i = m/(m + mo), and fiQ = 1 — [i = mo/(m + mo). The temperature ratio 7 is determined from 
the condition 



_ S r Q,xy 

1 f* P* ' 

SO r xy 

where the "cooling rate" Q — Co/ 11 f° r the impurity is given by 



d-l 



1/2 



(l + ao) [l-|(l + 0)(l + a o ) 



(A16) 



(A17) 



For elastic collisions {a = a = 1 ), Eqs. (TA"3)) - (|A"5|) and (|A10p - (|A15p lead to P£. = P *^ = S ir In this case, as 
expected, the solution to Eq. (|A16[) is 7 = 1. In addition, if we assume that particles of the gas and the impurity are 
mechanically equivalent (i.e., m = mo, a — 00, and a — ao), then P*j = PJy , C* = Co an d so 7 = 1. Beyond these 
two limit cases, the temperatures T and Tq are different so that there is a violation of energy equipartition. 
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APPENDIX B: BEHAVIOR OF THE ZEROTH-ORDER PRESSURE TENSORS NEAR THE STEADY 

STATE 



This Appendix addresses the behavior of the pressure tensors P*j and Pq i • of the gas particles and the impurity, 
respectively, near the steady state. The behavior of the second-degree moment P*j was studied in Ref. [2fJ in the 
three-dimensional case (d = 3). The extension to an arbitrary number of dimensions is straightforward and here only 
the final expressions are displayed. The behavior of the yy-element is given by 



dPl 



4P* 



<A + P* y 



da* J s yy 2afA + d(2f3 + (* 
where a* (a) is the steady state value of a* (p, T) given by Eq. (|A6|) , 



(Bl) 



1 + a 



1 - 



1 



2d 



-{I -a) 



and A = (dP*/da*^ is the real root of the cubic equation 

d 2 

2a* s i A 3 + Adaf(C + /3)A 2 + — {7(* 2 + 14C0 + 4/3 2 )A + d 2 /3((* + P)~ 2 {W 2 - X* 2 ~ PC)- 

In the above equations, it is understood that all the quantities are computed in the steady state. 
Let us consider now the elements of Pq^- In dimensionless form, they verify the equation 



1 1 * 9 



p* 



a ilPo,jl 



B* 



where 



B* 



m 
x Q vp 



|d CCiCj - J [/ (0) ,/ ( °)]. 



(B2) 



(B3) 



(B4) 



(B5) 



Upon deriving Eq. (|B4|) , use has been made of the fact that in the hydrodynamic regime the dimensionless pressure 
tensor depends on p and T only through its dependence on the reduced shear rate a* — a/v(p,T) [see Eqs. (|55[) and 
([56]) ] . The collisional moment B*j can be estimated by using Grad's approximation (1A8|) with the result PJ| 



B* — YSij + XqPq^j + XP* 



where 
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Let us consider the elements Pq xy and P n * . From Eq. (|B4[) . one gets 

2 a da~* 



^ a Pxy + C* ) ( 1 _ 9° ) Po,xy + a *P(},yy ~ X Q PQ xy ~ XP xy , 



(BIO) 



(Bll) 
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- faPZv + C*) (l - P lyy - XoP * yy = Y + XP; y . (B12) 

This set of equations can be written as 

dPp,xy _ r, (C* + 2 a *Pxy + Xo)Pp,xy ~ a *Pp,yy + XP* y 



da* " a* (C* + §a*P 5 



(B13) 



xy) 



^o.yy _ 2 ^* + i a * P 4 + Xo ) P o,m + Y + m 14 x 
3a* a* (C* + 2a*P*J 

It must be remarked that the temperature ratio 7 = To /T is also a function of a* in the hydrodynamic solution. Since 
Pq,xx + [d— 1)Pq yy = dy, the corresponding equation for the derivative (dy/da*) can be obtained from Eq. (IB4j) as 

2 So- <> + |o»PJ s 

where £q is given by Eq. (j A17|) . The set of three coupled equations (|B13I) . (|B14[) . and (|B15|) has a singular point 
corresponding to the steady state solution, i.e., when a*(p,T) = a* (a). In this limit (a* — *■ a*), the numerators and 
denominators of Eqs. (|B13p . (|B 14|) . and (|B15|) vanish. The limit can be evaluated by means of l'Hopital's rule. Thus, 
when one differentiates with respect to a* the numerators and denominators of Eqs. (|B13|) . (|B14|) . and (|B15[) and 
then takes the limit a — ► a* , one gets the relations 

/ OP* \ /1 \ _1 

' " 1 [^-alXs ~ *o] {XsPo.xy + (KPlxy + Ky) (dj/da*) s 

- [Plyy + a* s (dP^ yy /da*) s ] + X(dP* y /da*) s } , (B16) 
(<>'\.,,\ Xs Po, vy + {V + X' G Pl yy + X'P; y ) (d>y/da*)s + X(0P; y /0a*) s 



(B18) 



V da* 



V da* ) s \a* sXs -X 

dy\ _ Xsl ~ § [Po,x y + a*JdP^ xy /da*) s 



where the subscript s means that the derivatives are computed in the steady state. In addition, we have introduced 
the quantities 

x s = \ + a *s( dP * v / da *)s] > ( B19 ) 

and Y' = (dY/dy), X' = (OX/dj), X' Q = (OY/dj), and Co = (dQ/dy). Note that all the quantities appearing on the 
right-hand side of Eqs. (|B16|) . (|B17j) . and (|B18j) are evaluated in the steady state. This set of equations can be easily 
solved to give the corresponding derivatives. The result is 



where 



d2_\ = Ai 
da* ) A 2 ' 



1 . „ \ f/1 . „ \ ( 2 

2 



(B20) 



A i = d [ 2 a *sXs -X Q \U -a* sXs -X \ I Xsl - j P o,xy 



< [XsP^xy ~ Plyy + X (0 P* y / '0 ^) ,] | 

+2af [xsPl yy + X(OP; y /da*) s ] , (B21) 

A 2 = d (~a* sX s ~ X \ (±a* sXs ~ *o) ((* + \xs + 7 C ) 

+ ~a* s (X' P^ xy + X'P* yy ) - 2a? (Y' + X' P$ m + X'P*) . (B22) 
Substitution of Eq. (|B20p into Eqs. (|B16|1 and (|B17|1 gives the derivatives {d a *P^ xy )s and (9 a »P * yy ) s , respectively. 
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APPENDIX C: CHAPMAN-ENSKOG -LIKE EXPANSION 



In this Appendix, we provide some technical details on the determination of the first-order approximation by 
means of the Chapman-Enskog-like expansion. Inserting the expansions (f2"9"| and ([51]) into Eq. (f2"5)) , one gets the 
kinetic equation for /, 



(i) 



^ + (V + u s ).V + ^.A 
mo ov 



f 



(o) 



(CI) 



The velocity dependence on the right side of Eq. (IC1[) can be obtained from the macroscopic balance equations 
l^nj) to first order in the gradients. They are given by 



d^xo = -(u s + 5u) ■ Vx , 



(C2) 



d { t 1} Su = -(u s + <Su) • VJu - - (v ■ P (0) - «f) , 



(C3) 



d[ 1] p = -(u s + 5u) • Vp - p(V • Su + C (1) ) - ^ (aPW + P (0) : Vdu) , 



(C4) 



d t (1) T = -(u s + Su) • VT- A (aPW + p( 0) : V«5u) ~K (1) , 
where p = toti is the mass density, 

4 1} = / <fcroc iC ,/«(c), 

and 



= — \- I <lc 1 1 ic- 
dp 



(j[/ (0) ,/ (1) ] + J[/ (1) ,/ (0) ]) 



Use of Eqs. ((C2I) dC5j) in Eq. ([CT]) yields 



t ] - aV v ^\ / (1) - JoI/o^./^^Ao-Vrro + Bo-Vp 

+C • VT + D : V5u + E • ^ 

+C (1) (^ + T9 r )/r + J [/ (0) ) /W], 



(C5) 



(C6) 



(C7) 



(C8) 



where we have introduced the force J- = Fo — (mo/m)F. The coefficients of the field gradients on the right-hand side 
of (|C8j) are functions of c and the hydrodynamic fields. They are given by 



g f (0) 
OXq 



(C9) 



g A0) ldf (°)dP {0) 

i>0,i(Cj — C; H 



(CIO) 



C 0li (c) 



s/f 1 ac dpi 

dT Cl p dSuj dT ' 



(Cll) 



(0) 



<3p 



dSiii 



dp dT J 



(C12) 
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J_dfl 

m da 



(0) 



(C13) 



Upon writing Eq. (|(J12|) use has been made of the expression of the pressure tensor P^. of the gas (2 j 



(1) _ _ dSu k 

Uj 'Hjkt g > 



(C14) 



where r]ijk£ is the viscosity tensor. Moreover, the expression for the distribution derived in Ref. (20| is given by 
Eq. OB- 

The solution to Eq. (|C8[) has the form 



(i) 



(C15) 



The coefficients Aq, Bq, Cq, T>q, and £q are functions of the peculiar velocity and the hydrodynamic fields xq, p, T, 
and Su. The cooling rate depends on space through its dependence on p and T. Moreover, there are contributions 
from (9j ' acting on the pressure, temperature, and velocity gradients given by 



'MP® , 2C (oA Vp _f^^_l^ VTi 



d dp 



d dT 2 T 



(C16) 



a (0) vr = 



v ('f „p£> + c«» 



dp 2 
2a 
dp 



1 + T 



d_ 

df 



\ p(0) i I^(o) 
/ xy 2^ 



Vp 
VT, 



(C17) 



(C18) 



The corresponding integral equations (|42| - (|46|) can be obtained when one identifies coefficients of independent gradi- 
ents in (|C8P and takes into account Eqs. (|C16|) - (|C18|) and the mathematical property 



dl 0) X = 



dX 
dp 



91 P 



OX 
dT 



a (0) T- 



dX 
dSui 



di 0) 6u t 



(C19) 



where in the last step we have taken into account that X depends on Su through c = V — 5u. 
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